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In this paper, we investigate the asymptotic nature of the quasinormal modes for 
"dirty" black holes — generic static and spherically symmetric spacetimes for which 
a central black hole is surrounded by arbitrary "matter" fields. We demonstrate that, 
to the leading asymptotic order, the [imaginary] spacing between modes is precisely 
equal to the surface gravity, independent of the specifics of the black hole system. 

Our analytical method is based on locating the complex poles in the first Born 
approximation for the scattering amplitude. We first verify that our formalism agrees, 
asymptotically, with previous studies on the Schwarzschild black hole. The analysis is 
then generalized to more exotic black hole geometries. We also extend considerations 
to spacetimes with two horizons and briefly discuss the degenerate-horizon scenario. 
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1. Introduction 

It has long been known that small perturbations of a black hole spacetime will 
evolve according to a one-dimensional Schro dinger- like equation [1]. Although no 
(normalizable) bound-state solutions are known, there is still a very interesting class 
of solutions, commonly referred to as quasinormal modes [2, 3, 4, 5, 6], that can be used 
as a "basis" for investigating the physics. The premise behind the quasinormal modes 
is to impose "radiation boundary conditions" at the edges of the spacetime (typically, 
spatial infinity and the black hole horizon); a constraint which necessitates a complex 
value for the frequency (and, hence, the "quasi" nomenclature) [7]. As it turns out, 
this complex frequency can be labeled by a discrete quantum number (n = 0, 1,2, ... ) 
and takes on the following asymptotic form [8]: 

kqnm{n) —in (gap) + (offset) + 0[n~ l l 2 ] as n — > oo , (1) 

where the "gap" and "offset" are model-dependent, complex parameters that are 
determined by the precise form of the "potential barrier" in the Schrodinger equation. 

For illustrative purposes, let us suppose an asymptotically flat spacetime, in which 
case the gap is purely real. A useful physical interpretation is that, for a highly damped 
black hole (as described by asymptotically large n), the real part of the offset measures 
the frequency of emitted radiation, whereas the gap corresponds to quantized increments 
in the inverse relaxation time [7] . Hence, one might expect that both of these quantities 
should be fundamental properties of a given black hole. In this sense, it is interesting to 
consider the modes for the specific case of a scalar or gravitational perturbation outside 
of a Schwarzschild black hole (as generated numerically by Nollert [9], Andersson [10], 
and substantiated analytically by Motl and Neitzke [11, 12, 13]): 

k^(n)=i^(n+\) + ^ + 0[n-^ as oo , (2) 

where m is the black hole mass and G is Newton's constant. | For future reference, let us 
re-express this result directly in terms of the surface gravity at the horizon (k = l/[AGm] 
for a Schwarzschild black hole): 

kqnmin) = in (^n + ^ + k + 0[n~ l/2 } as n — > oo . (3) 

Note that the gap in the Schwarzschild case is precisely the surface gravity. One 
might be reasonably inclined to wonder if this occurrence is an artifact of a particularly 
simple model (having only one dimensional parameter) or, rather, a resilient feature 
of more exotic black hole geometries. Undoubtedly, this is an important question but, 
perhaps strangely, there has been limited progress towards finding the answer. Studies 
on quasinormal modes have, for the most part, been of a highly model-specific nature 
(either numerical [9, 10, 14, 15] or analytical [11, 12, 13, 16, 17]). 

| Here and throughout, all other fundamental constants are set to unity and a 3+1-dimensional 
spacetime is presumed. 
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In the current paper, we directly address the above question with a generic 
analytical approach that is elegant in both its simplicity and its general applicability. 
(This continues a long-term project of one of the current authors if black 
hole thermodynamics is as fundamental as expected, then it should be completely 
independent of the detailed specifics of the black hole under consideration. The 
thermodynamic features of the black hole system should be generic and unaffected 
by any "dirt" that might surround the black hole [18, 19, 20].) In our analysis, we 
utilize the first Born approximation for the scattering amplitude [21] as a means for 
approximating the reflection coefficient. It is then possible to identify the quasinormal 
modes by locating the poles in the scattering amplitude. (See, for example, [13].) Note 
that the Born approximation is normally thought of as a "high-energy" approximation. 
More specifically, it becomes increasingly accurate as \k\ — > oo (this is undoubtedly 
true for real k and, in all likelihood, complex k as well). Although this approach may 
appear "quick and dirty", it does certainly confirm the known asymptotic behavior in 
the Schwarzschild case. Moreover, the basic methodology can readily be extended to 
generic black hole spacetimes. We ultimately show that the gap does indeed equate with 
the surface gravity; this being a model-independent feature of black hole geometries. 

Before proceeding, let us point out that there has been, quite recently, substantial 
discussion within the gravity community concerning the real part of the quasinormal 
frequency (i.e., the real part of the "offset") [22, 23, 24]. This has followed, to a large 
degree, from Hod's observation [25] that the offset can be used to fix, uniquely, the 
spacing between the black hole area eigenvalues. (This observation implies that the 
area of a black hole horizon can be quantized to give an evenly spaced spectrum; a 
notion that was first advocated by Bekenstein [26].) Unfortunately, we are not able 
to address this sub-leading (offset) term at the level of the first Born approximation. 
Nonetheless, if quasinormal modes do indeed play some sort of role in quantum gravity 
(in spite of the fact that they arise out of a purely classical process), then it becomes an 
important matter to not only calculate but physically motivate the complete spectral 
form. We would argue that the findings of this paper — specifically, identifying a 
universal property of quasinormal mode spectra — can be viewed as progress in this 
direction. 

The rest of the paper is organized in the following manner. The next section 
concentrates on the quasinormal modes of the Schwarzschild black hole, as this 
(relatively) simple case nicely illustrates the formalism and allows for a direct comparison 
with known results. In Section 3, we go on to consider the quasinormal modes for 
"generic" black hole spacetimes (the only stipulations being staticity and spherical 
symmetry), including those with more than one horizon. At the end, we also discuss 
the conceptual difficulties of a degenerate-horizon scenario. Section 4 provides a brief 
summary and discussion. 
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2. Schwarzschild black holes 



Although our formal techniques have quite general applicability, the focus will be, 
for the time being, on the minimally complicated (but nonetheless interesting) case 
of a Schwarzschild black hole. We are, simplistically speaking, interested in small 
perturbations of the spacetime outside of the black hole horizon. It is therefore 
appropriate to begin with a standard Klein-Gordon equation, 

-^d„ gT d v ] tt(r, t,9,4>)=Q, (4) 



V-9 

where \& is the massless perturbation field and the metric describes the background 
(four-dimensional Schwarzschild) spacetime. That is, 



,hr I l —) ( lf 2 ■ ( 1 



2m V 



, dr 2 + r 2 dn 2 . (5) 
r J \ r J 

(Note that, for the remainder, we set G — 1. If need be, G and other fundamental 
constants can be easily re- introduced via dimensional considerations.) 
Employing a separation-of-variables technique and writing § 



#(r,0,<M) = - Hr) Y tm {d,<j>) eM+ikt) , 
r 

one can convert equation (4) into a Schro dinger- like equation of the form [1] 
d 2 

— ip - V[r(n)] ip = -k 2 i) , 

where is the so-called tortoise coordinate defined by [27] 
dr* _ 1 
dr 1 — 2m/r 

For scalar perturbations, the "scattering potential" is found to be 

2m\ 



V(r) = (l - ^) 



+ 1) 2m 



(6) 



(7) 



(8) 



(9) 



where t is the orbital angular momentum (/ = 0, 1, 2, ...) and let us emphasize that V(r) 
is a rational polynomial in r but not in r*. There is a natural generalization to higher 
spin fields; namely [28] , 



V(r) 



1 - 



2m 



+ 1) , 2m(l-j 2 ) 



(10) 



and, for the physically most relevant cases, 



1 

3/4 


-5/4 
-3 



scalar 
Dirac 

vector 



Rarita-Schwinger j 
gravity j 



3 = 
3 = 1/2 

3 = 1 
3/2 

2 . 



(11) 



§ With the assumption that the imaginary part of k is positive — cf, equation (2) — we can fix the 
sign in the temporal exponent by requiring the modes to be exponentially decaying in time. 
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r*(r) = r + 2m In 



(12) 



In more generic circumstances, the (generalized) tortoise coordinate will be a 
complicated and perhaps unsolvable function of r. In the present case, however, 
r* = r*(r) can readily be obtained from a straightforward integration of equation (8). 
This process yields 

"r — 2m" 

. 2m . 

Although this function is not explicitly invertible, we do know that the region r G 
(2m, oo) maps into the region r* e (— oo, +oo). That is to say, the exterior of the black 
hole maps into the entire real line, and we effectively have a one- dimensional scattering 
problem. 

It is natural, at this point, to impose the physical boundary conditions of purely 
outgoing plane waves at spatial infinity (r* — > oo) and purely ingoing plane waves at the 
horizon (r* — > — oo). Such boundary conditions can only be realized for complex values 
of k, thus leading to the notion of quasinormal-mode solutions to the wave equation [2]. 
(For the physical relevance of these modes, see the previous section; see also [7].) 

For a generic (one-dimensional) scattering problem, it is known that the 
quasinormal modes can be identified with the poles of the reflection coefficient or, 
equivalently, the scattering amplitude as a function of complex asymptotic momenta [7]. 
Here, we will approximate the scattering amplitude by way of the first Born 
approximation [21]. Generally speaking, this approximation is obtained from the Fourier 
transform of the scattering potential with respect to the momentum transfer. More 
specifically, we can regard 

/+oo 
V[r(r*)} exp[+2iAr»] dr* (13) 
-oo 

as an approximate form of the scattering amplitude. Take note of the factor of +2 in the 
exponential. This follows from the fact that the momentum transfer is always minus 
twice the incident momentum in a one-dimensional scattering: 

q = k f -k i = (-k i )-k i = -2k i . (14) 

Our objective is now clear: identify the poles in equation (13) with V(r) as given 
by equation (9). As an initial step, let us perform a change of variables, 

r+°° dr* 
a(k) oc / V[r]exp[2ikr*(r)] — — dr , (15) 
J 2m dr 

and then incorporate equations (8) and (12) to give 

\r - 2m] i4mk ' 1 
2m 

A trivial shift in the integration variable then leads to 

r+oo 2m + r r 1 *4mfc— 1 



a(k) oc / V[r] - — exp[2i/cr] 

J2m 2m 



dr . (16) 



r+oo Am + r 

a(k) oc / y [2m + r] — exp[4ikm] exp[2ikr] 

Jo 2m 

For future convenience, let us now define 



2m 



dr . (17) 



z = -2ikr that is r = iz/(2k) , (18) 
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and so obtain 



a(k) oc — 



-i iAmk 



X 



I 



Amk 

+00 



V 



exp [4ikm] 

(Amk + iz) exp[-z] z^' 1 dz. 



(19) 



2m+ 2k 



At this point, one can, with a brief inspection of the integrand, anticipate the presence 
of a linear combination of Gamma functions. || For instance, a (hypothetical) constant 
term in V(z) would necessitate contributions that go as T (iAmk) and T (iAmk + 1). The 
identification of poles then becomes a trivial exercise; in this case, k qnm (n) = in/ '(4m) , 
where n is any non-negative integer. 

Let us now be more precise and recall our specific form (9) for the scattering 
potential. It follows that 



V 



2m + Tk 



IZ 



Amk + iz 



t(i + 1) 2m(l - j 



2m+f kJ 



V 



+ 



(2m + 



iz \ ' 
2k J 



Substituting into equation (19), we then have 

iAmk 

exp [4ifcm] 



a(k) oc - 



x 



Amk 

+00 



(20) 



(21) 







1} 



2m(l — j 2 ) 



exp[-z] z iAmk dz . 



( 2m + *i) ( 2m + ^ k j 

Since the only immediate concern is the location of the poles, it is sufficient, for our 
purposes, to regard z as a small parameter and Taylor expand the quantity in the square 
brackets. Such an expansion yields a power series in ascending powers of z, and it soon 
becomes evident that we obtain 



a(k) oc ]T C s T(i4mk + 1 + s) , 



(22) 



s=0 



where each of the coefficients, C s , is itself a regular and well-defined quantity which we 
do not need to calculate. 

We can now readily locate the poles in the scattering amplitude. They occur at 

iAmk + 1 = — n where n > , (23) 

so that, up to the validity of the first Born approximation (which presumably gets better 
as I A: I becomes larger), we have 

n > . (24) 



kqnm(jl) 



n 



i- — where 
4m 



Let us re-emphasize that the first Born approximation can only be expected to have 
validity for very large scattering energies; meaning that this result can only be trusted 

|| Keep in mind the primary definition of the Gamma function for positive real n: T(n) = 
Jo°° s " _1 ex P I~ s ] ds- Its analytic continuation to the complex plane has poles at all non-positive 
integers. 
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when n » 1 (or really when n — > oo). Fortunately, this is just the regime we are 
interested in.^f 

Our outcome for the location of the quasinormal modes agrees, asymptotically, 
with the results obtained from both numerical studies [9, 10] and other analytical 
means [11, 12]. That is, 

Ti 

k qnm (n) = i- hO[l] when n — > oo. (25) 

This spacing between asymptotic modes, or the gap, can alternatively be written as 
l/4m = k, where k is the surface gravity of the Schwarzschild black hole. At a first 
glance, this could be interpreted as just a happy coincidence; insofar as there is only one 
dimensional parameter in the problem. (That is, the gap would almost certainly have 
to be the surface gravity times a numerical factor.) It will, however, be shown below 
that the asymptotic spacing is universally given by the relevant k, irrespective of the 
details of the black hole spacetime. 

Before proceeding, let us point out that this result for the Schwarzschild gap has 
been known, even analytically, for quite some time. For instance, Liu and Mashoon 
[29], as well as Andersson [30], have made this observation by identifying the high- 
frequency limit of the scattering equation with the confluent hypergeometric equation. 
It is, however, unclear how one would translate this method, or any of the more recent 
analytical treatments [11, 12], into a generic setting. 

3. Generic black hole spacetimes 

In this section, the previous formalism will be extended to "generic" black hole 
spacetimes. By generic, we mean static and spherically symmetric but with an otherwise 
arbitrary geometry. (So the black hole can be "dirty" in that it may be surrounded by 
an arbitrary source of static and spherically symmetric matter [18, 19, 20].) Eventually, 
we will elaborate on spacetimes with two (or possibly more) horizons, with the second 
horizon not necessarily being an event horizon per se (e.g., the popular model of a 
Schwarzschild black hole enclosed by a de Sitter cosmological horizon, the Kottler or 
Schwarzschild-de Sitter geometry). + Horizons are (for the time being) assumed to be 
non-degenerate; but the special case of horizon degeneracy (e.g., an extremal Reissner- 
Nordstrom black hole) will be discussed at the very end of the section. 

f Indeed the Born series converges in relatively few cases, and is more typically an asymptotic series. 
Therefore we only expect the Born approximation to be sensitive to the leading order contributions, at 
least when applied naively. 

+ Although the formalism can also be extended to multiple-horizon scenarios, it is the feeling of the 
authors that, from an operative viewpoint, a single observer in any spherically symmetric geometry 
would be able to deduce the existence of at most two horizons. That is, for any such observer, the 
accessible spacetime will have no more than two spatial boundaries and, therefore, no more than two 
(non-degenerate) horizons can ever come into play. Hence, formal considerations will be restricted to 
spacetimes with one or two horizons. 
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3.1. Single- Horizon Scenarios 



Given a static spacetime and spherical symmetry, the metric for a generic black hole 
can always be expressed, without loss of generality, in the following manner (see, for 
example, [18, 31]): 

-l 



ds 2 



-24>{r) 



2m(r) 



dt 2 +[l 



2m(r) 



dr 2 + rW. 



(26) 



Here <p{r) is a model-dependent function (related to the Morris-Thorne "redshift 
function" [32]), while the "mass" parameter m(r) is equivalent to the Morris-Thorne 
"shape function" [32]. Alternatively, we can write 



ds 2 



-2<P(r) 



2m(r) 



-dt 2 + 



ar 



2m(r) 



r 2 dQ 2 



e -2*(r) (i 

which leads, quite naturally, to a generalized tortoise coordinate, 
dr* 1 



and so 



dr e -Hr) (i _ 2 ™( r )) ' 
2m(r)' 



ds 2 = e~ 2 ^ 



1 - 



-dt 2 + dr 2 +r 2 dQ 



(27) 



(28) 



(29) 



Next, let us next revisit the Klein-Gordon equation (4) for a massless perturbation 
field. In terms of the generic spacetime described above, this equation expands into 

2m(r)\ 



\-d 2 t +y 2 d rJ 2 d rf +e- 2 ^) (l- 



# = , 



(30) 



where A2 represents the angular part of the d'Alembertian. We can now proceed in 
the standard way; namely, factoring \1/ into a temporal part [exp(+ikt)], angular part 
(the usual spherical harmonic) and a radial part. One can express the radial part as 
r p u[r(r*)]; in which case, the unique choice of p which eliminates all terms containing 
d rt u (but not the double derivatives) happens to be p = —1. (This was exactly the 
same exponent as occurred in the Schwarzschild case, which is now seen to be generic.) 
Given this choice, a Schrodinger-like form is once again obtained: 

^ u -V[r(n)]u = -k\ (31) 
dr; 

except that the scattering potential is significantly more complicated than was found 
for the Schwarzschild case. More specifically, some straightforward calculation yields 
the following result:* 



V(r) = - (d 2 j) + e 



-2<«r) 



2m(r) 



+ 1) 



-2<p(r) 1 _ 



2m(r) 



+ 1) 



~f(r) 



(32) 



(33) 



* Here, for sake of simplicity, we are considering a scalar (j = 0) perturbation. It is, however, technically 
possible to extend the calculation to arbitrary j, and such a generalization should not alter any of our 
results or conclusions. 
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where we have defined 



f(r) = 1 



2m(r)\ <j>'(r) | 2 / _ m(r)' 



(34) 



and a prime indicates a derivative with respect to r. 

Provided that /(r) is a well-defined and regular quantity (which must always be 
the case since this is equivalent to the tortoise coordinate being well behaved), it is 
clear that the generic potential is qualitatively very similar to that of the Schwarzschild 
scenario [cf, equation (9)]. Rather, the sticking point with the current calculation is 
that there is, in general, no means of obtaining a closed form solution for r* = r*(r). 

We can, however, circumvent the forementioned difficulty by first taking note of 
the exact generic expression for the surface gravity [18], 



K 



1 

2r~ h 



-4>(r h ) 



1 - 8np(r h ) r\ 



1 _d_ 

2 d7 



-<t>(r) 



1 - 



2m(r) 



(35) 



Th 



where r = r\ indicates the horizon [specified by 2m{r\) = rh] and p{r) is the energy 
density. Hence we can expand equation (28) to give 
dr* 1 
dr 



2« {(r - r h ) - a(r - r h ) 2 + 0[(r - r h f]} 
!1 ' r h)] , 



Yk + ^ 



(36) 



2K{r - r h ) 

where a is a model-dependent constant parameter. (Although such an expansion is 
technically valid only in the vicinity of the horizon, it turns out to be sufficient for 
the purpose of identifying poles in the scattering amplitude. As seen in the previous 
section, the location of the poles, at the level of the first Born approximation, is only 
sensitive to the near- horizon geometry.) Moreover, the above expansion can now be 
directly integrated to yield 



a , 
2k 



1 , 
H In 

2k 



r-r h 



+ 0[{r-r h ) 



(37) 



Beginning with the (approximate) expression for the scattering amplitude (13) and 
repeating the steps that took us up to equation (19), we now find that (here using the 
"convenient" choice of z = —iakr/K and neglecting the irrelevant prefactors) 



a(k) oc 



V 



IZK 

rh^ -r 

ak 



- + iz + 0[z 2 } ) exp[-z] (l + 0[z 2 }) z l(kjK) - x dz 



(38) 



Now what about the scattering potential? An inspection of equation (33) reveals that 
V [rh + izK/(ak)] will translate into a power series in z times several terms; each of 
which contains some (typically negative) power of the argument [77, + izn/(ak)] . As 
in the previous section, we can, for our purposes, regard z as small and Taylor expand 
appropriately. Hence, we end up with a complicated power series in z. Fortunately, the 
precise details of this expression are unimportant to us. 

Given the previous expression (38) and the above discussion, it becomes clear that 
r +00 

(39) 



a(k) oc / exp[-z] z tk/K F(z) dz , 
Jo 
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where F(z) is some unknown but, in principle, calculable power series in integer powers 
of z. [Given that the surface gravity is non-zero, the leading term in F is of the order 
z°, as can be seen by evaluating dV(r)/dr ^ at the horizon.] 
From the above form, it follows that 

00 ( k \ 

a(k) «J]C ( r i- + l + s , (40) 

s =o V K J 
(As before, when it comes to locating the position of the poles, we do not need to 
calculate the coefficients C s .) Hence, the first-Born-approximated poles are located at 

k qnm {n) = inn where n > , (41) 

which is enough to imply that the actual physical poles lie at 

kqnm(n) = inn + 0[l\ where n > , (42) 

with the constant term becoming (in comparison) irrelevant as n — > 00. This is the 
main result of the paper and substantiates our claim at the end of Section 2. 

3.2. Dual-Horizon Scenarios 

Let us now envision an observer "trapped" between two horizons; for instance, the region 
of spacetime between a Schwarzschild black hole and a de Sitter cosmological horizon 
or the two event horizons in a Reissner-Nordstrom black hole (although our formalism 
certainly allows for much more generic situations). We propose that it is straightforward 
to extend the previous analysis to such situations by virtue of the following observations: 

(i) Given a spherically symmetric and static spacetime, the metric can still be cast, 
without loss of generality, in the form of equation (26). 

(ii) In the vicinity of any given horizon, the derivative dr\j dr (and all related quantities) 
can be expanded in the form of equation (36); where the surface gravity (k) and 
horizon location (r h ) are uniquely defined parameters for the horizon in question. 

(iii) At the level of the first Born approximation, the poles in the scattering amplitude 
are sensitive only to the near-horizon geometry (or geometries) of the spacetime. 

To further elaborate, starting with the appropriately revised form of equation (15), 

f T2 dr* 
a(k) oc / V[r] exp[2i£x*(r)] — — dr (43) 
Jn dr 

(where T\ and r 2 locate the two horizons in the spacetime), we can split this integral 

at some intermediate point (say r x ) and then, by way of the third observation, make 

the approximation of extending the integrals to r x — > ±00. This procedure effectively 

yields two distinct sets of poles, one coming from each horizon. (These are actually 

two distinct scattering problems which can be distinguished by the orientation of the 

incident wave.) That is, one can anticipate the asymptotic form (with the surface 

gravities labeled accordingly) 

kqnminx) = % n,\ k x + 0[l\ or k qnm (n 2 ) = % u 2 K 2 + 0[l] , (44) 
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as nior2 - ► oo. Note that the inverse of the surface gravity effectively fixes the time 
scale, so that the modes scattered by the inner/outer horizon (i.e., the larger /smaller 
surface gravity) will dominate observations at earlier/later times. Further note that 
this phenomena agrees with analytical estimates, based on the use of the Poschl-Teller 
potential, as performed by Suneeta [33]. 

Implicit in the above discussion is that the two horizons can indeed be spatially 
isolated. This immediately rules out spacetimes with degenerate (or very nearly 
degenerate) horizons; for instance, an extremal Reissner-Nordstrom black hole or a 
Nariai (degenerate Schwarzschild-de Sitter) spacetime. To further complicate matters, 
degenerate horizons invariably have a vanishing surface gravity, thus rendering our 
previous expansions to be useless in this context. 



3.3. Extremal horizons 



In spite of the difficulties inherent to degenerate-horizon spacetimes, one possible 
recourse would be to focus on the near-horizon form of the second radial derivative 
of g tt (since the first derivative vanishes at a degenerate horizon). That is to say, we can 
now expand dr*/dr as follows: 



dr* 
dr 



1 



_a(r - r h ) 2 + fi(r - r h ) 3 + ^(r - r h ) 4 + 0[(r - r h ) 5 ]_ 
1 (3 , f3 2 -ia 



+ 



+ 0[(r-r h )} , 



(45) 



a(r — rh) 2 a 2 (r — rh) o.° 

where a (related to the second derivative), (5 and 7 are model-dependent constants. The 
tortoise coordinate is then 



a[r — rh) a z 



r -r h 



+ 



(r-r h )+0[(r-r h ) 2 } (46) 



Now, closely following the previous methodology, we find that [after shifting the r 
integration from (r h , 00) to (0, 00)] the scattering amplitude takes on the form 

2ik 1 I3 2 
+ — 

a r 



a(k) oc 



V [r h + r] exp 

r ^ -2ik/3/a 2 

r h 



7a 



2ikr + 0\r' 



x 



1 pi 



+ 0[1] \ dr 



(47) 



With the definition z = r/r^, the above can be rearranged into 



a(k) oc 



+00 



exp 



2ik 1 P 2 ~ja n ., 
1 5 — 2ikr h z 



ciTh z a° 

where F(z) is again some (in principle) calculable power series in z with integer 
exponents. 

An evaluation of the z integral now yields modified Bessel functions which, apart 
from the possibility of a trivial pole at k — 0, do not have poles at finite \k\. We can 
physically interpret this finding as follows: When the surface gravity is non-zero, the 
poles in complex /c-space can ultimately be traced back to the small-z (or near-horizon) 
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behavior of the integrand. Whereas, in the current case of a degenerate horizon, the 
integrand factor exp(— ik/r) washes out any possibility of a pole for k ^ 0. 

Which is to say, by setting k — ■» in the general analysis, we are able to "predict" 
that the quasinormal modes all collapse to zero momentum. But what this really means 
is that the quasinormal modes (if any exist) do not lie in a region of the complex k 
plane where the first Born approximation is trustworthy. Thus, what we really expect 
for extremal black holes is that the quasinormal modes (if any exist) lie in some bounded 
region of the complex k plane and are either finite in number or densely scattered in 
some bounded region. 

Before concluding, let us point out some discrepancies between our analysis and 
other recent works that have considered the quasinormal modes of a charged black 
hole; both analytically [13, 16] and numerically [14]. Firstly, there is evidence that 
the gap of a Reissner-Nordstrom black hole — if it is periodic at all — goes, not as 
the surface gravity at the outer horizon, but rather as a complicated function of both 
surface gravities (inner and outer horizon). Naively, this is not what one would expect 
from our findings, given that the interior horizon is not explicitly part of the scattering 
problem. Secondly, the same papers have found that the highly damped modes of an 
extremal Reissner-Nordstrom black hole are formally identical to a Schwarzschild black 
hole of the same mass. This is clearly contrary to our arguments above. It is quite 
possible that our methodology breaks down, in some subtle way, in spacetimes with a 
"hidden horizon"; that is, it is feasible that the first Born approximation is unable to 
properly account for these added complexities. It should be noted, however, that some 
of these authors have also commented on interpretative difficulties [16] and numerical 
instabilities [14] on account of some "peculiar features" [34] of the quasinormal spectra 
of charged black holes. Moreover, the Schwarzschild limit of their Reissner-Nordstrom 
spectrum does not appear to give back the desired result. We would suggest that the 
quasinormal mode problem for this charged model requires further attention. 

4. Conclusion 

To summarize, we have used a simple analytical method — based on the first Born 
approximation for the scattering amplitude — to locate the quasinormal modes for 
a static, spherically symmetric but otherwise generic black hole spacetime. For the 
very special case of a Schwarzschild black hole, our methodology was found to agree, 
asymptotically, with the results obtained from various numerical and analytical studies. 
Moreover, that the asymptotic spacing or "gap" is equivalent to the surface gravity (as 
is accepted in the Schwarzschild case) was shown to be a model-independent feature of 
a wide range of black hole spacetimes. We also generalized this outcome to spacetimes 
with two horizons and reported some limited progress in the context of degenerate- 
horizon geometries. 

The virtues of our formal treatment include generality and a straightforward 
analytical approach that does not obscure the physical process being investigated. 



Dirty black holes: Quasinormal modes 



13 



Nevertheless, the simplicity of our method has a price: we have (so far) only been 
able to locate the leading term in the asymptotic expansion of the quasinormal mode. 
Meanwhile, the next-order term — that is, the "offset" - - has sparked considerable 
recent interest because of a conjectural relationship with the black hole area spectrum. 
It is unfortunate that the Born approximation seems to diverge at higher orders, so that 
our approach would have to be significantly modified to reproduce the entire spectrum. 
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